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ABSTRACT 



5_j ' Aims. Results from y-ray observations by the H.E.S.S. telescope array in the direction of the young stellar cluster Westerlund 2 are presented. 
I Metiiods. Stereoscopic imaging of Cherenkov light emission of y-ray induced showers in the atmosphere is used to study the celestial region 
around the massive Wolf-Rayet (WR) binary WR 20a. Spectral and positional analysis is performed using standard event reconstruction tech- 
niques and parameter cuts. 

Results. The detection of a new y-ray source is reported from H.E.S.S. observations in 2006. HESS J1023-575 is found to be coincident 
with the young stellar cluster Westerlund 2 in the well-known HII complex RCW 49. The source is detected with a statistical significance of 
more than 9a-, and shows extension beyond a point-like object within the H.E.S.S. point-spread function. The differential y-ray spectrum of the 
emission region is measured over approximately two orders of magnitude in flux. 

Conclusions. The spatial coincidence between HESS J1023-575 and the young open cluster Westerlund 2, hosting e.g. the massive WR binary 
WR 20a, requires one to look into a variety of potential models to account for the observed very-high-energy (VHE) y-ray emission. Considered 
emission scenarios include emission from the colliding wind zone of WR 20a, collective stellar winds from the extraordinary ensemble of hot 
and massive stars in the stellar cluster Westerlund 2, diffusive shock acceleration in the wind-blown bubble itself, and supersonic winds breaking 
out into the interstellar medium (ISM). The observed source extension argues against a single star origin of the observed VHE emission. 

Key words. ISM: star forming regions - ISM: individual objects: HESS JI023-575, WR 20a, RCW 49 (NGC 3247, G284.3-0.3) - gamma- 
rays: observations 
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1. Introduction 

The prominent giant HII region RCW 49 (also NGC 3247), 
and its ionizing cluster Westerlund 2 are located towards 

t now at School of Chemistry & Physics, University of Adelaide, 
Adelaide 5005, Australia 
* now at DESY Zeuthen 
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the outer edge of the Carina arm of our Milky Way. 
RCW 49 is a luminous, massive star formation region, 
and has been extensively studied at various wavelengths. 
Af ter the initial report based on ph otographic images 
bv lRodgers. Campbell & Whiteoak (Il960l) . the region was ob 



serve d at X- ray energ ies with Einstei n ([Hertz & G: 



1984 



Goldwurm. Caraveo & Bignamil 



indlav 
1987h . 



ROSAT dBelloni & Mereghettilll994 . and recently Chandra, 
which discov ered ~500 po i nt sou rces in the vicinity of 



RCW 49 jTsuiimoto et al.1 I2004D 
spatially coinci dent with the 
Westerlund 2 (iTownsley et al.l l2005l) 



with -100 of them 
central open stellar cluster 
Mid-infrared mea- 



surements with Spitzer revealed still ongoin g massive star 



formation in RCW 49 dWhitnev et al.1 120041) . The regions 



surrounding Westerlund 2 appear evacuated by stellar winds 
and radiation. The surrounding dust is distributed in fine 
filaments, knots, pillars, bubbles, a nd bow shocks through 



out the rest of the HI I complex (iChurchwell et al.l 12004 : 
Conti & Crowthen 120041) . Radio continuum observations by 



ATCA at 1.38 and 2.38 GH z indicate two wind-blown shells 
in the core of RCW 49 dWhiteoak & Uchidal Il997l) : one 
surrounding the central region of Westerlund 2, the other the 
prominent WR star WR 20b. The distance to RCW 49 is still 
uncer tain and values r ange between ~2.2 kpc ferand & Blitzl 
1993b up to 7.9 kpc dMoffat. Shara & Potteilll99lh . whereas 
intermediate values o f 4.2 kpc from 21 cm absor ption line 
profile measurements ( McClure-Griffiths et alj|2001 ). 5.75 kpc 
from the distance estimate to wards the prominent WR star 
WR 20a dVan der Huchtll200lh. and 6.4 k pc from photometric 



measurements ( Carraro & Munari 20041) corresp ond to more 
recent refinements. Finally, iRauw et alj d2007l) presented a 



compelling re-determination of the distance to Westerlund 2 
by spectro-photometric measurements of 12 cluster member 
O-type stars of (8.3+ 1.6) kpc. This value is in very good 



agreem ent with the (8.0+1.0) kpc as measured by iRauw et al 



d2005l) . determined from the light curve of the eclipsing 
binary WR 20a. In summary, we adopt the distance value 



of the weighted mean of (8.0+1.4) kpc dRauw et al.1 l2007n 



throughout the manuscript, thereby associating WR 20a as a 
cluster member o f Westerlund 2. 

Westerlund 2 dWesterlund 1 19601) contains an extraordinary 
ensemble of hot and massive stars, presumably at least a dozen 
early-type O-stars, and two remarkable WR stars. Two Wolf- 
Rayet (WR) stars have been found within RCW 49, WR 20b 
at a2ooo = 10*'24™18.H, dim) = -57°48'30", and WR 20a 
at g?_ nnn = 10''23"58.'^0, ^2000 = -57°45' 49" dVander Hucht 
2001b . both first reported as WR stars by IShara etalj dl99lk . 
the latter believed to be a member of the stellar cluster 
Westerlund 2. Only recently w as WR20a estabUshed to be a 
binary: both lRauw et alJ d2004l) and lBonanos et al.l d2004l) pre- 
sented solutions for a circular orbit with a period of 3.675, 
and 3.686 days, respectively. Based on the the orbital period, 
the minimum masses have been found to be (83 + 5) M© and 
(82 ± 5)M p for the primary and secondary components, re- 



spectively dRauw et al.ll2005h . This classifies the WR binary 



WR 20a as the most massive of all confidently measured bi- 
nary systems in our Galaxy. An orbit al inclination of ~75 ° was 
inferred from photospheric eclipses ( Bonanos et alj|2004 ). and 



the spectral classification was s ubsequently refine d as a binary 
consisting of two WN6ha stars dRauw et al.ll2005l) . The super- 
sonic stellar winds of both WR stars collide, and a wind-wind 
interaction zone forms at the stagnation point with a reverse and 
forward shock. In a detached binary system like WR 20a, the 
colliding wind zone lies between the two stars, and is heavily 
skewed by Coriolis forces. The wind velocity at the stagnation 
point is likely lower than i n other positions due to th e radia- 
tive inhibition mechanism dStephens & Pollocklll994l) caused 
by the companion's stellar photon field. In a close binary sys- 
tem the winds can only be radiatively accelerated to a fraction 
of their expected Voo ~ 2800km/s, and a co mparatively low pre - 
shock wind velocity of ~ 500km/s follows dRauw et al ] l2005h . 

WR stars are commonly thought to represent an evolved 
phase for the most massive stars in our Galaxy, however the 
most luminous hydrogen-rich WR stars m ay still be thought o f 
as the most massive, main-sequence stars ( Moffat et al. 20061) . 
exhibiting strong, broad WR-like emission lines at optical 
wavelengths. Some WR stars are already established as non- 
thermal radio emitters, whereas claimed associations with y- 
ray emission remain to be confirmed. It has been suggested 
that non-thermal radio emission observed in a significant frac- 
tion of the massive stars could be relate d to colliding winds 
in binaries ( Dougherty & Williams 2000l) . An important ques- 
tion is therefore whether all massive stars which produce 
non-thermal radio emission are indeed binaries. Another rele- 
vant issue is whether the relativistic electrons which produce 
the observed non-thermal (synchrotron) radio emission also 
produce detectable (inverse Compton) high-energy emission. 
Synchrotron emission has not been detected from the WR 20a 
system, presumably because of free-free-absorption in the op- 
tically thick stellar winds along the li ne of sight. However, 
WR 20a has been detected in X-rays ( Mereghetti & Bellonil 
19961) . With a flux of the order of 10 erg cm^^ 5-\ the 
soft X-ray emission is consistent with the predicted inten- 
sity for the thermal radi ation from the shock-h eated mate- 
rial in the stellar wind dPittard & StevensI 120021) . However, 
non-thermal and thermal components of the X-ray emission 
are currently indistinguishable. In high-energy y-rays, three 
unidentified EGRET sources have been found in the wider 
vicinity of RCW 49, 3EG J1013-5915 , 3EG J1014-5705, and 
3EG J1027-5817 dHartman et al.ll999l) . These EGRET sources 
are characterized by both source confusion, and indication of 
extended emission, therefore viable counterparts are not estab- 
lished readily. 

It has been suggested that in massive binary systems 
high-energy y-rays are produced either by optically-thin 
inverse Compton scattering of relativistic electrons with 
the dense photospheri c stellar radiation fields in the wind 



wind c ollision zone ( Eichler & Usov| 1993 *: Whit e & Cher 



1995 



Benaglia & Romero 



20061: IPittard & Doughertvl I2TO 



iReimer. Pohl & Reimei 
or in hadronic emis- 
sion scenarios. Here, y-rays are either directly produced 
in neutral pion decays, with the mesons produced by 
inelastic interac tions of rela t ivistic nucleons with the 
wind material dWhite & ChenI II992I: iBenaglia & Romero 



Torres et alj I2OO 



20031: 

Pittard & Doughertvl |200€ 



IReimer. Pohl & Reimeil 12006 : 
or alternatively are produced 
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in cascade models ( IBednarekl l2005h . In the latter scenario, 
y-radiation originates from inverse-Compton pair cascades, 
which are initiated by high-energy neutral pion decay photons 
(from nucleon-nucleon interactions in the stellar winds). Due 
to the dense radiation field at the wind colUsion zone in the 
close binary WR 20a, y-rays of a few tens of GeV are optically 
thick to photon-photon absorption, and pair cascades are 
unavoidable. If those nuclei reach sufficiently high energies 
(despite the expected weak shock in the extremely dense 
radiative environment), the projectile nucleons themselves 
are produced by photodisintegration of he avy nuclei in th e 
radiation fields of the massive WR-stars ( Bednarekl 12005 ). 
For a high magnetic field of order lO^G, particles could 
be energized up to ~ lO'^eV (adopting typical WR-wind 
parameters) in the magnetic field reconnection scenario, and 
up to 10^ GeV if first-order Fermi acceleration takes place 
there. The maximum power extractable from WR 20a must be 
a fraction of the kinetic wind energy provided by the system, 
estimated to be ~ 10^'^erg/s. 

Detectable VHE y-radiation from the WR 20a binary sys- 
tem is so far only predicted in the pair cascade model, although 
detailed modeling of the WR 20a system in the other scenar- 
ios is still pending. For VHE y-rays, photon-photon absorption 
will diminish the observable flux from a close binary system 
such as WR 20a ( lBednarekl2 005; Dubus 2006). 

Here we report on the first VHE y-ray observations towards 
the WR binary WR 20a, and its immediate surroundings with 
the H.E.S.S. telescopes. The observations and analysis yield 
the detection of a new source, HESS J1023-575. This source 
is found to coincide with the central region of RCW 49, and 
accordingly, with the stellar cluster Westerlund 2, which hosts 
the WR binary WR 20a. The results and possible interpretation 
of these findings are presented in the following. 

2. H.E.S.S. Observations and Data Analysis 

The H.E.S.S. (High Energy Stereoscopic System) collabo- 
ration operates a telescope system of four imaging atmo- 
spheric Cherenkov telescopes located in Namibia (23° 16' 17" 
S 16°29'58" E), at 1800m above sea level. Each telescope has 
a tesselated spherical mirror with 13 m diameter and 107 m- 
area, and a telescope spacing of 120 m. Cherenkov light emit- 
ted by extended air showers in the atmosphere is imaged by 
a high resolution camera (960 photo-multipliers, pixel size 
0.16°). Shower parameters and the pr imary particle type are 



2004 



determ ined by image reconstruction dAharonian et al. 
2006bh . The telescopes are operated in coincidence mode. 



wobble-mode observations, where the telescopes are pointed 
offset from the nominal source location to allow for simultane- 
ous background estimation. The wobble offsets for these obser- 
vations range from 0.5° to 2°, with the majority of data taken 
with wobble offset less than 0.8°. The zenith angles range be- 
tween 36° and 53°, resulting in an energy threshold of 380 GeV 
for the analysis. 



-57° 




which requ ires a trigger fro m at least two out of the four 



10^30'" 10''25'" 10''20'" 

Right ascension (J2000) 

Fig. 1. Correlated excess sky map with an oversampling radius 
of 0.12° of the Westerlund 2/RCW 49 region. The background 
in each bin is estimated using a ring around the test position. 
The map has been smoothed with a two-dimensional Gaussian 
of radius 0.04° to reduce the effect of statistical fluctuations. 
Overlaid contours correspond to statistical significances of 5, 
7, and 9cr, respectively, determined using the nominal oversam- 
pling radius for extended sources of 0.2°. The inlay in the lower 
left corner shows the excess distribution for a point-source de- 
rived from Monte Carlo data with the same zenith-angle and 
offset distribution as the data. The cross denotes the best fit 
position of the gamma-ray data with Icr statistical uncertain- 
ties. The WR binary WR 20a is indicated by an upright filled 
triangle in the Westerlund 2 stellar cluster (dashed circle), the 
reversed fiUed triangle denotes the location of WR 20b. 

The data have been analyzed using the H.E.S.S. standard 
Hillas analysis with standard cuts (image amplitude > 80 p. e.). 
O ptical efliciency correctio ns have been applied as described 



telescopes (IFunk et al.l 120041) . The H.E.S.S. telescope array in lAharonian et al.l (l2006bl) . A point source analysis on the 



achieves a point source sensitivity above 1 TeV of < 2.0 x 



10" cm" s" (1% of the flux from the Crab nebula) for a 5cr 
detection in a 25 hour observation. 

The dataset described here consists of 14h (12. 9h live 
time) of data taken between March and July 2006, either on the 
nominal source location of WR 20a or overlapping data from 
the ongoing Galactic plane survey. Quality selections were im- 
posed on the data, excluding those taken during bad weather 
or with hardware irregularities. The data have been obtained in 



nominal position of WR 20a resulted in a clear signal with a 
significance of 6.8 standard deviations. Further investigations 
of the skymap of this region revealed an extended excess with 
a peak significance exceeding 9 standard deviations (Fig. 1). 
The position of the center of the excess was derived by fitting 
the two-dimensional point spread function (PSF) of the instru- 
ment folded with a Gaussian to the uncorrected excess map: 
ffiooo = 10'>23™18^ + 12\ ^2000 = -57°45'50" + 1'30". The sys- 
tematic error in the source location is 20" in both coordinates. 
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Fig. 2. Number of excess events versus the squared angular dis- 
tance from the best fit position of the excess {0-). The dashed 
Une shows the expectation for a point source derived from 
Monte Carlo data. The solid line is a fit of the PSF folded with 
a Gaussian (cr = 0.18° + 0.02°). 




E ( TeV ) 



Fig. 3. Differential energy spectrum and fit residuals to a sin- 
gle power-law fit of HESS J1023-575 from photons inside the 
85% containment radius (0.39°) around the best fit position. 
The background is estimated with background regions of the 
same size and distance from the camera center as the signal 
region. The spectrum can be described by a power law with a 
photon index of F = 2.53 +0.16stat±0.1syst and a normalization 
at 1 TeV Oo = (4.50+0.56stat+0.90syst)xlO-'"TeV-i cm-^g-'. 
The integral flux for the whole excess is F(E>380GeV) = (1.3 
+ 0.3) xl0-"cm-2s->. 



The number of excess events versus the squared angular dis- 
tance to this position in comparison to the expectation for a 
point source (dashed line) is shown in Fig. 2. The excess is 
clearly extended beyond the nominal extension of the PSF. A 
fit of a Gaussian folded with the PSF of the H.E.S.S. instru- 



ments gives an extension of o" = 0.18° + 0.02°. The differen- 
tial energy spectrum for photons inside the corresponding 85% 
containment radius of 0.39° is shown in Fig. 3. It can be de- 
scribed by a power law (dN/dE= <l>o • (E/1 TeV)"'^) with a pho- 
ton index of F = 2.53 + 0.16stat ± O.lsyst and a normalization at 
1 TeVof Oo = (4.50+0.56stat+0.90syst)xlO-i-TeV-i cm-^g-'. 
The integral flux for the whole excess above the energy thresh- 
old of 380 GeV is (1.3 + 0.3) xlO"" cm ^ g-i. No significant 
flux variability could be detected in the data set. The fit of a 
constant function to the lightcurve binned in data segments of 
28 minutes has a chance probability of 0.14. The results were 
checked with independent analyses and were found to be in 
good agreement. 

3. HESS J1 023-575 in the context of y-ray 
emission scenarios 

The detection of VHE photons from the region 
RCW 49/Westerlund 2, characterized by a moderately 
hard power-law spectrum with an index F ~2.5 is indicative of 
the existence of extreme high-ene rgy pa r ticle acceleration pro - 
cesses in this young (~2-3 Myrs: IPiattil ( Bica & Claria 1998 )) 
star forming region. The estimated luminosity above 380 GeV 
is ~ 1.5 X lO^^erg/s at a distance of 8 kpc, which corresponds 
to 0.5% of the total kinetic energy available from the colliding 
winds of WR 20a, and ~1.6% of the kinetic energy of WR 20b. 
With a projected angular size of submilliarcsecond scale, the 
WR 20a binary system, including its colliding wind zone, 
would appear as a point source for observations with the 
H.E.S.S. telescope array. Unless there are extreme diff'erences 
in the spatial extent of the particle distributions producing 
radio. X-ray, and VHE y-ray emission, scenarios based on 
the colliding stellar winds in the WR 20a binary system face 
the severe problem of accounting for a source extension of 
0. 18° in the VHE waveband. At a nominal distance of WR 20a 
(8.0kpc; Rauw et al. 2007), this source extension is equivalent 
to a diameter of 28 pc for the emission region, consistent 
in size with theoretical pre dictions of bubbles blown from 
massive stars into the ISM dCastor et al. 1975 ). The spatial 
extension found for HESS J1023-575 contradicts emission 
scenarios where the bulk of the y-rays are produced close to the 
massive stars. If VHE y-rays were indeed produced near the 
massive stars in the binary, a regular m odulation of th e y-ray 
flux due to absorption would be expected lDubusI (l2006h . Future 
observations of HESS J1023-575 will allow one to probe for 
such a pattern, which could discard or unambiguously label 
the origin of the VHE emission. 

Alternatively, the emission could arise from collective ef- 
fects of stellar winds in the Westerlund 2 cluster. Diffusive 
shock acceleration in cases where energet ic particles experi- 
ence multiple shocks ( Klepach et al. 2000l) can be considered 
for Westerlund 2. The stellar winds may provide a sufficiently 
dense target for high-energy particles (accelerated to hundreds 
of TeV in the winds themselves or driven by supernova explo- 
sions), allowing the productio n of Tr^-decay y-ray s via inelastic 
pp-interactions. According to lRauw et al.l (120071) Westerlund 2 
may host about 4500 Mq in the form of stars with M > Mq. 
The mechanical energy injected through the steUar winds de- 
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termines mainly the total energy available for collective wind 
effects, estimated to be about 5 .4 x IQ-^^erg/s at most. A collec- 
tive w ind scenario proposed bv lDomingo-Santamaria & Torres 
( 20061) predicts that the extension of a y-ray source corre- 
sponds to the volume filled by the hot, shocked stellar winds. 
Variability is generally not expected since flux modulation of 
the winds of individual stars or binaries is believed to average 
out in a collective wind. Interestingly, a putative VHE source 
may not possess a non-thermal counterpart at MeV/GeV en- 
ergies in such a scenario, since convection will prevent low- 
energy particles from entering the wind. In both colliding wind 
zone and collective wind models, one would expect sub-TeV 
spectra that resemble those of Supernova remnants. 

Magneto-hydrodynamic (MHD) particle acceleration (e.g. 
by multiple shocks or turbulence produced by supersonic 
flows) in a magnetized plasma may also be considered when 
particles penetrate into a dense medium. In such a case, the 
massive stellar winds of Westerlund 2 could accomplish three 
tasks: they could ensure sufficient particle injection into the tur- 
bulent plasma, feed the magnetic turbulence with energy via 
wind-wind interactions in the massive star association, or pro- 
vide copious photons and dense material to serve as target for 
TeV-photon producing particle-photon and photon-photon in- 
teractions. Supershells, molecular clouds, and inhomogeneities 



dBykov & Toptyginlll987l) . Shocks and MHD turbulent motion 
inside a stellar bubble or superbubble can efficiently transfer 
energy to cosmic rays if the particle acceleration time inside the 
hot bubble is much shorter than the bubble's expansion time. 
X-ray observations of RCW 49 revealed, apart from point-like 
emission from individual stars, a diffuse component centered 
on the ionizing cluster, widely ext ending outside the core, with 
a soft component and a hard tail (iBelloni & Mereghettilll994 ; 



Townslev et al.ll2005l) . This indicates the existence of a wind- 
blown shell from the stars at the cluster core, filled with a 
hot (~ 0.1-3 keV) tenuous plasma shocked by turbulence. 
Two wind-blown bubbles are revealed from continuum radio 
observations: one around WR 20b with a diameter of ~ 4.1' 



and an expansion velocity of 67 km/s (IShara et al.llI991l) . the 
other aroun d the Westerlund 2 cluster center with a diameter 
of ~ 7.3' jWhiteoak & Uchidal [l997l) . A ridge of enhanced 
GHz radio emission lies between these two shells, and it has 
been suggested that th is may represent a site of collision be- 
tween these two shells ( Whiteoak & Uchidalll997l). Most in ter- 
estingly, the extent of a blister dWhiteoak & Uchidal 1997 ) on 
the western side of the Westerlund 2 bubble appears to be com- 
patible in direction and location with the center of gravity of 
HESS J1023-575. Such a blister is indicative for rapid expan- 
sion into a low-density medium outside the wind-blown bubble. 
Figs. 4 and 5 show overlays of the detection significance con- 
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Fig. 4. HESS J1023-575 significance contours (corresponding 
5, 7 and 9 cr), overlaid on a B-band image from the Second 
Palomar Observatory Sky Survey (POSS-2). The filled circle 
denotes the best fit position with Icr statistical uncertainties. 
The WR binary WR 20a is indicated by an upright filled tri- 
angle in the Westerlund 2 stellar cluster (dashed circle), the 
reversed filled triangle denotes the location of WR 20b. 

embedded in the dense hot medium may serve as the targets 
for y-ray production in Cosmic Ray interactions. Such en- 
vironments have been already studied in the nonlinear the- 
ory of particle acceleration by large-scale MHD turbulence 
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Fig. 5. HESS J1023-575 significance contours (corresponding 
5, 7 and 9 cr), overlaid on a 843 MHz image from the Molonglo 
Observatory Synthesis Telescope (MOST). Black symbols and 
dashed circle as in Fig. 1. The wind-blown bubble around 
WR 20a, and the blister to the west of it are seen as depres- 
sions in the rad io continuum map. The bli ster is indicated by 
white dots as in I Whiteoak & Uchida ( 1997 ). and appears to be 
compatible in direction and location with the center of gravity 
of HESS J1023-575. 



tours of HESS J1023-575 (E>380 GeV) on the B-band image 
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of the region as obtained from the Second Palomar Observatory 
Sky Survey (POSS-2), and on the 843 MHz radio continuum 
image from the Molonglo Observatory Synthesis Telescope 
(MOST), respectively. Whereas the open cluster Westerlund 2 
does not perfectly coincide with the peak of the emission seen 
at very high energy y-rays (Fig. 4), the direction of the outbreak 
and spatial location of the blister is a better match (Fig. 5). The 
later association also relaxes constraints arising of the detection 
of HESS J1023-575 as an extended VHE y-ray emitter. 

Shock acceleration at the boundaries of the blister may en- 
able particles to diffusively re-enter into the dense medium, 
thereby interacting in hadronic collisions and producing y- 
rays. Consequently, a scenario as outlined in Volk (1983)) for 
a supernova-driven expansion of particles into a low density 
medium may be applicable to the expanding stellar winds into 
the ambient medium. If one accepts such a scenario here, it 
might give the first observational support of y-ray emission 
due to diffusive shock acceleration from supersonic winds in a 
wind-blown bubble around WR 20a, or the ensemble of hot and 
massive OB stars from a superbubble in W esterlund 2, breakin 
out beyond the edge of a molecular cloud /Tenario- Tagle'l97 
Volk & Forman 1982: Cesai-sky & Montmerle .1983; Bykov i 



200lh . Accordingly, one has to consider that such acceleration 
sites will al so cont ribute to the observed flux of cosmic rays in 
our Galaxy dCasse & Paul 1980). 

Other open star clusters hosting WR-stars are known to 
supersede the total injected mechanical energy from winds 
of massive stars as available in Westerlund 2: Westerlund 1, 
which is older than W esterlund 2 and hes at a recently re- 
vised distance of 5 kpc dCrowflier et al.ll2006l) . and NGC 3603, 
a very prominent galactic starburst regio n at a distance of 
7 kpc. Both offer approxim ately 10^''erg/s ( Muno et al. 2006t 



EP 

i 



Stephens & Hartwelll 120031) . Search for TeV-emission from 



open stellar clusters has been pre viously carried out for a list 
of northern clusters by HEGRA ( lAharonian et alj|2006cl) . re- 



sulting in upper limits in the order of 10"^' to lO"'"' cm"^ s"' at 
varying energy thresholds between 0.8 and 3 TeV The south- 
ern hemisphere stellar cluster Westerlund 2 was not targeted 
by HEGRA. The detection of very-high-energy y-rays associ- 
ated with Westerlund 2 by H.E.S.S. clearly motivates system- 
atic searches for TeV-emission from massive open stellar clus- 
ters. 



4. Conclusions 

H.E.S.S. observations have led to the discovery of an extended 
source HESS J1023-575 in the direction of Westerlund 2. 
Possible sources are the massive WR binary system WR 20a 
(although source extension and variability studies do not sup- 
port a colliding wind scenario at present), the young stellar 
cluster Westerlund 2 (although the cluster itself appears too 
compact to account for the observed VHE emission in a col- 
lective wind scenario), and cosmic rays accelerated in bubbles 
or at their termination shock and interacting with their environ- 
ment. Further observations with the H.E.S.S. telescope array 
will help to discriminate among the alternatives in the interpre- 
tation of the observed VHE y-ray emission. 



However, the convincing association of HESS J1023-575 
with a new type of astronomical object — a massive HII re- 
gion and its ionizing young stellar cluster — profoundly distin- 
guishes this new detection by the H.E.S.S. telescope array al- 
ready from other source findings made during earher G alactic 
Plane Scan observations ( Aharonian et al ] l2005ll2006ah . 
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